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We analyze the scaling behavior of the higher Lyapunov exponents at the Anderson transition. 
We estimate the critical exponent and verify its universality and that of the critical conductance 
distribution for box, Gaussian and Lorentzian distributions of the random potential. 



I. INTRODUCTION 

Universality is a concept that is central to the the- 
ory of critical phenomena and continuous phase transi- 
tions. The idea is that critical phenomena are indepen- 
dent of the "details" of the physical system depending 
only on its symmetry and dimensionality. All systems 
in the same "universality class" should exhibit the same 
critical phenomena!!! 

The Anderson transition is a continuous phase transi- 
tion and so the concept of universality should apply. In 
this case, the important symmetries are those of time re- 
versal and spin rotation. Three universality classes are 
expected: orthogonal, unitary, and symplectic. In this 
paper we are concerned with the orthogonal universal- 
ity class, which is comprised of systems that have both 
time reversal and spin rotation symmetry. This corre- 
sponds to an electron in a random potential in the ab- 
sence of a jamgnetic field and with negligible spin-orbit 
inter actionntl 

One of the simplest ways to test the validity of the con- 
cept of universality is to change the distribution of the 
random potential in which the electron moves. Since this 
does not change the universality class, the critical expo- 
nents for different distributions of the random potential 
should be the samea. In a recent papera we presented 
numerical estimates of the critical exponent for the An- 
derson transition for three different distributions of the 
random potential: box, Gaussian and Lorentzian. The 
estimates were made by analyzing the finite size scaling 
of the localization length of electrons on long wires, and 
the results support the contention that the critical expo- 
nent has a universal value independent of the distribution 
of the random potential. 

In this paper, wie analyze the scaling of the higher Lya- 
punov exponentscl for the box, Gaussian and Lorentzian 
distributions of the random potential. We give a 
clear demonstration of single parameter scaling for these 
higher Lyapunovs, a clear verification of the universality 
of this scaling with respect to these three distributions 
of random potential, and show that this scaling yields 
estimates of the critical exponent that are in almost ex- 
act agreement with our previous estimates based on the 
scaling of only the lowest Lyapunov exponents. (The low- 



est Lyapunov exponent is the inverse of the localization 
length.) 

We also present an estimation of the critical con- 
ductance distribution p c (g)- This is expected to 
take a size independent universal form at the critical 
pointQ. Thepsize independence of p c (g) has been verified 
previouslyHii3. Here, we verify that the form of p c (g) is 
independent of the distribution of the random potential 
for the three distributions mentioned. 



II. SCALING OF THE HIGHER LYAPUNOV 
EXPONENTS 

A. Model and method 

The motion of the electrons is described by the Ander- 
son model 



(1) 



<%,j> 



where Cj(Ci) denotes the creation (annihilation) opera- 
tor of an electron at the site i of a three dimensional cubic 
lattice. The value of the potential at site i is Wj . Hop- 
ping is restricted to nearest neighbors and its amplitude 
is taken as the unit of energy, V = 1. 

We have examined three different distributions of the 
random potential. 

1. A box distribution with each Wi uniformly dis- 
tributed on the interval [-W/2,W/2\. 

2. A Gaussian distribution with mean zero and vari- 
ance W 2 /12. 

3. A Lorentzian distribution of the form 

W 



p(Wi) 



(2) 



for which the variance diverges. 



The standard method used i& transfer matrix stud- 
ies of the Anderson transition^ involves rewriting the 
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Schroedinger equation for electrons on a long bar with 
dimensions L x L x L z as a product of transfer matrices 



M = JJ M n , 

n=l 

where M„ is the transfer matrix, 
a matrix Q by 

n = In mm\ 



(3) 



From M we can define 



(4) 



with eigenvalues {— i/jvj • • • , — vi, ■ ■ ■ , +^jv} where 
N = L 2 and the i/s are arranged in increasing order. 
The eigenvalues occur in pairs of pnppsite sign as a con- 
sequence of current conservationli3tL3 and therefore we 
need only consider the positive v' s. The Lyapunov ex- 
ponents {a n } for the random matrix product Q3J) can be 
defined as 



a n = lim = — 



(5) 



These were calculated to within a specified accuracy us- 
ing the method described inllil. The localization length A 
for electrons on the bar is the reciprocal of the smallest 
positive Lyapunov exponent A = \/ot\. In our previ- 
ous work, we examined only the scaling of the quantity 
A = X/L = \ju\L. Here we have examined the scaling 
of the quantities z n defined by 



a n L, 



(6) 



for higher values of n. In performing the scaling analysis 
corrections to scaling were taken into account using the 
method described in (Ref. 5). We assumed that the data 
obey the following scaling law 



(7) 



where tf> is a relevant scaling variable and <f> is an irrele- 
vant scaling variable. We approximated this scaling func- 
tion by its first order expansion in the irrelevant scaling 
variable and fitted the data for each z n to the form 



(8) 



The scaling variables were approximated by expansions 
in terms of the dimensionless disorder w = (W c — W) /W c 
where W c is the critical disorder separating the insulating 
and metallic phases, 



(9) 



The critical exponent v describes the divergence of the lo- 
calization length (correlation length on the metallic side) 
as the transition is approached 



i = Co M 



(10) 



Note that the absolute scale of the localization length 
£o cannot be determined with the method we have used 



here. The decay of the irrelevant scaling variable with 
system size is described by an exponent y < 0. 

The functions F (n) and F[ n) were expanded either to 
second or third order in w. Also, when possible, the 
expansion for tp was truncated at the first order in w. 
The decision at which order to truncate these various se- 
ries was based on the criteria of obtaining an acceptable 
goodness of fit with as few parameters as possible. 



B. Results 

We simulated systems with sizes L = 
4,5,6,8,10,12,14 and 16 and Fermi energy E F = 
calculating data to 0.1% accuracy with the exception 
of some data very close to the critical point that were 
calculated to 0.05% accuracy. We found that the z„'s 
obey single parameter scaling provided that the system 
sizes considered are large enough. The minimum system 
size required increased with n so that, for example, we 
could fit only the data with L > 10 when analyzing z\q 
for the box distribution. In practice, we could demon- 
strate scaling for the first ten exponents for the box and 
Gaussian distributions. For the Lorentzian distribution 
the required system sizes seem to be larger and we could 
demonstrate scaling only for the first three exponents. 

The res ults of the scaling analysis are presented in Ta- 
bles H - III and some of the data are plotted in Figs. |l|j^. 
The estimates of the critical disorder obtained by scaling 
different Lyapunov exponents for a given distribution of 
the random potential are in very close agreement. The 
same can also be said for the critical values of the pa- 
rameter z n for various n between different distributions 
of the random potential. The critical values z^ for a 
given n is defined by 



z(") = F o (n) (0), 



(11) 



and is expected to be universal if the scaling function 
for each Lyapunov is universal. The data clearly sup- 
port this idea. Turning to the critical exponent ^, we 
find estimates that are remarkably consistent across all 
the Lyapunov exponents and disorder distributions we 
considered. 



III. THE ZERO TEMPERATURE 
CONDUCTANCE DISTRIBUTION 

Next we considered the zero temperature conductance 
distribution when W = W c , i.e. the critical conductance 
distribution at the Anderson transition. We set the Fermi 
energy Ep = and, based on the results of Section II, 
we used W c = 16.54 for the box distribution, W c = 21.3 
for the Gaussian distribution, and W c = 4.26 for the 
Lorentzian distribution when simulating the conductance 
distribution. An ensemble of samples of cubic shape was 
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generated, and for each sample, a Green's function iter- 
ation technique was used to determine the tpuismission 
matrix t in a two probe measuring geometryO-3. The di- 
mensionless conductance g was calculated using the Lan- 
dauer formula 

g = 2trtt\ (12) 

where the factor of two takes account of spin degeneracy. 
(We would like to emphasize that the t matrix that ap- 
pears in this formula is an N x N matrix, not an L 2 x L 2 
matrix, where N is the number of propagating modes at 
the Fermi energy in_the contacts attached at the left and 
right of the samplePj) . 

In Fig. U and |(], the conductance distribution for en- 
sembles of 12 x 12 x 12 cubic samples for the three distri- 
butions of random potential are presented. Fixed bound- 
ary conditions were imposed on the wavefunction in the 
direction transverse to current flow so that the system 
size L = 12 is large enough that the distribution shown 
are a good approximation to the true critical distribu- 
tions PcU&) and p c (\ng) obtained in the asymptotic limit 
L — > ooc3. (Larger system sizes are required if periodic 
boundary conditions are considered i-and, the asymptotic 
form of the distribution is different£Elll3'EZl. ) Note that 
despite the appearance of Fig. [|, p c (g) — ► as g — > 0. 
This can be deduced from the form of p c ( lna) .m Fig. ^ 
or a more detailed plot of p c (g) near g = OEStLZl. 

Looking at Figs. [| and || we see that critical distri- 
butions obtained for the three different distributions of 
the random potential are almost identical. This is con- 
sistent with the assertion that the critical conductance 
distribution is universal. 



IV. SUMMARY AND DISCUSSION 

The critical exponent for the Anderson model with 
random bopping but no random potential has been 
estimated^ to be v = 1.61 ± .07. Since models with 
random hopping and random potentials are in the same 
universality class, we should, and do, obtain estimates for 
the critical exponent that are consistent with Ref. 18. 

The computer time required to estimate the Lyapunov 
exponents with the method used heretil scales as L 7 for 
a given accuracy. Thus, it is unlikely that this method 
can be extended to much larger systems sizes in the near 
future. The critical exponent of the Anderson transition 
has also been estimated by-analyzing the scaling of en- 
ergy level statistics (ELS)Ej. The most recent, and in 
our opinion the most accurate, estimate of v obtained 
with ELS is v = 1.45 ± 0.2E9. This is consistent with 
our results but the estimated uncertainty is an order of 
magnitude worse than we have achieved here. 

In the simulations whose results we have presented 
here, the Fermi energy is always set at the band cen- 
ter Ep =0. For one dimensional systems it has been 
demonstrated that single parameter scaling docs not hold 



for states with energies in a band gap of the cocceaoond- 
ing ordered (i.e. zero disorder W = 0) systemETo . In 
that case, single parameter scaling is only observed for 
sufficiently strong disorder. We suspect that this might 
also hold for the higher Lyapunov exponents in the three 
dimensional systems we have studied here for energies 
outside the band. 

To conclude, we have demonstrated that the higher 
Lyapunov exponents obey single parameter scaling at 
the Anderson transition and presented evidence that 
this scaling behavior is universal and independent of the 
choice of distribution of the random potential. Analysis 
of this scaling behavior yields estimates of the critical 
exponent that are remarkably consistent between differ- 
ent Lyapunov exponents and different distributions of the 
random potential. Finally, we presented evidence that 
the critical conductance distribution is also universal and 
independent of the choice of the distribution of the ran- 
dom potential. 
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n 


N d 


N p 


W c 


(n) 
Zc 


V 


1 


336 


12 


16.54(53,55) 


1.737(733,742) 


1.56(55,58) 


2 


336 


12 


16.55(54,55) 


2.801(797,804) 


1.57(56,58) 


3 


292 


12 


16.55(55,56) 


3.596(591,602) 


1.57(56,58) 


4 


250 


10 


16.56(55,57) 


4.25(24,26) 


1.58(57,59) 


5 


250 


10 


16.57(56,59) 


4.83(82,84) 


1.58(57,59) 


10 


162 


10 


16.59(51,71) 


7.0(6.9,7.2) 


1.55(52,59) 



TABLE I. The results of the scaling analysis for the box distributed random potential. Here n indicates which Lyapunov 
exponent is being considered, N d the number of data available, N p the number of parameters used to fit the data, W c is the 



estimate of the critical disorder, z c is defined in ([ll]) and v is the estimated critical exponent. The figures in brackets are 
95% confidence intervals for the corresponding critical parameter estimated using the bootstrap method (Ref. 14). 



n 


N d 




N p 


Wc 


(») 

Z c 


V 


1 


200 




10 


21.30(28,31) 


1.737(734,742) 


1.58(56,61) 


2 


200 




9 


21.29(28,30) 


2.796(794,799) 


1.57(56,58) 


3 


175 




9 


21.31(30,32) 


3.593(588,597) 


1.58(57,59) 


4 


150 




9 


21.31(29,33) 


4.244(236,252) 


1.58(56,60) 


5 


125 




9 


21.29(27,33) 


4.80(79,82) 


1.57(55,60) 


10 


100 




9 


21.40(30,57) 


7.05(6.98,7.18) 


1.55(50,62) 



TABLE II. The results of the scaling analysis for the Gaussian distributed random potential. The notation is explained in 
the caption for Table 1. 



n 


N d 


N p 


Wc 


(n) 

Zc 


V 


1 


256 


10 


4.25(24,26) 


1.70(68,71) 


1.57(52,61) 


2 


192 


10 


4.24(23,26) 


2.71(64,75) 


1.58(53,65) 


3 


127 


10 


4.29 


3.62 


1.58 



TABLE III. The results of the scaling analysis for the Lorentzian distributed random potential. The notation is explained 
in the caption for Table 1. (We could not execute the bootstrap method for n = 3 so no confidence intervals are available.) 
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FIG. 1. vs disorder W for the box distributed random potential. 
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FIG. 2. The same data as in Figure |l| after corrections to scaling have been subtracted and plotted versus L/£ where £ is 
the localization (correlation) length. 
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FIG. 4. The same data as in Figure after corrections to scaling have been subtracted and plotted versus L/£. 
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FIG. 5. The distribution of the conductance for ensembles of 100,000 12 x 12 x 12 cubes for box (solid line), Gaussian (dotted 
line) and Lorentzian (dashed line) distributed random potentials. 
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